Coarsening of granular clusters: two types of scaling behaviors 
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We report on an experimental study of small cluster dynamics during the coarsening process 
in driven granular submonolayers of 120^im bronze particles. The techniques of electrostatic and 
vertical mechanical vibration were employed to excite the granular gas. We measure the scaling 
exponent for the evaporation of small clusters during coarsening. It was found that the surface area 
of small clusters S vs time t behaves as S ~ (to — t) 2//3 for lower frequencies and S ~ (to — i] for 
higher frequencies. We argue that the change in the scaling exponent is related to the transition 
from three dimensional (3D) to two dimensional (2D) character of motion in the granular gas. 

PACS numbers: 45.70.Mg,45.70.Qj 



One of the most fascinating phenomena observed in 
freely-cooling and driven granular gases is spontaneous 
separation into dilute and dense regions jl], ^] . This phe- 
nomenon stimulated great interest among theorists and 
experimentalists J^, ||, pi pi . The clustering phenomenon 
appears to not be limited to just the traditional granu- 
lar systems with short-range hard-core collisions. It was 
shown that some granular systems with competing long- 
range electromagnetic interactions and short-range colli- 
sions exhibit somewhat similar behavior [Q. In the case 
of mechanically vibrated granular matter where the en- 
ergy is injected in the system through local collision with 
the bottom plate, clustering takes place due to the in- 
crease of the dissipation rate with the increase of the 
density of granular gas j| |). In the case of electro- 
statically driven granular matter there is an additional 
effect of electrostatic screening which decreases the ex- 
ternal electrostatic force acting on particles inside dense 
clusters jjj. 

The phase diagram for these two seemingly different 
systems appears to be rather independent on the specifics 
of driving [| [?J . Particles are immobile for small values 
of the driving force (amplitude of vibration acceleration 
in the case of mechanical system and the amplitude of 
electric field for the electrical system) . When the driving 
force exceeds the first critical value, isolated particles be- 
gin to move. If the driving force value is increased above 
a second threshold value, the granular medium forms a 
nearly uniform granular gas. Between these values the 
uniform gas phase is unstable and a fraction of particles 
can form immobile clusters. The clusters exhibit coars- 
ening dynamics of Ostwald-ripening type: small clusters 
disappear and large clusters grow with time. 

Here we report experimental study of small cluster dy- 
namics during the coarsening process in driven granular 
submonolayers. We systematically compare two tech- 
niques of energizing the granular gas: mechanical and 
electrostatic. We measure the scaling exponent for the 
evaporation of small clusters during coarsening. We have 



found that in both systems the surface area of small clus- 
ters S vs time t behaves as S ~ (to ~ t) 2 / 3 for lower fre- 
quencies and S ~ (to — t) for higher frequencies. This 
transition in scaling behavior coincides with the transi- 
tion from 3D dynamics of the granular gas to 2D dy- 
namics as the driving frequency is increased for a fixed 
acceleration. 

The experimental setup for mechanical shaking is sim- 
ilar to that in [|[ p], ||, ||] . We performed experiments in 
a rigid circular container of 13.5 cm diameter vibrated 
vertically by an electromagnetic shaker. The spheri- 
cal bronze particles (120 — 150/im) constituted approxi- 
mately 2/3 of a monolayer coverage on the optically-flat 
horizontal plate (single-crystal Silicon wafer). The ap- 
plied frequency varied in the range of 20-100 Hz. Most 
of the measurements were performed under atmospheric 
pressure in air. Some results were verified in vacuum. 

To excite a granular media electrostatically, particles 
were placed between plates of a large capacitor and en- 
ergized either by DC or AC electric fields Q. We used 
27 x 27 cm transparent capacitor plates (glass coated by 
indium doped tin oxide) and the plate spacing was 1.5 
mm. The particles used for these series of experiments 
are the same as for the mechanical shaking measurements 
and constituted approximately 1/3 of a monolayer cov- 
erage on the bottom plate. The applied field was varied 
in the range — 7 kV/cm and its frequency was varied 
from to 120 Hz. The experiments were performed in an 
atmosphere of dry nitrogen to reduce adhesion of parti- 
cles on the plate due to humidity in air Jh}. The CCD 
camera was suspended above both cells. We analyzed the 
total number of clusters, evolution of individual clusters, 
mean cluster sizes, etc. 

In the electro-cell, energy injection works as follows 
0. Particles in the contact with bottom plate acquire 
an electric charge and the electric force exerted on the 
particle is i*o = 1.36i? 2 £? 2 , R is the radius of the parti- 
cle, E is the external field. When Fo exceeds the grav- 
ity, the particle moves upwards, recharges upon colliding 
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with the upper plate and falls down. In the case of DC 
fields the process repeats in a cyclical manner. By ap- 
plying AC field one may turn the particle back before 
it reaches upper plate. Thus, increasing the frequency 
of the field on controls the extent of vertical motion of 
the particles. Clustering of granular gas with consequent 
ripening occurs if E\ < E < E% (E\ and E-% is the first 
and second threshold values). The clustering comes from 
two effects: increased dissipation of the kinetic energy in 
the regions with high granular media density and elec- 
trostatic screening (two close particles acquire smaller 
charge then two well-separated ones). In the mechanical 
cell, vibration amplitude plays the role of electric field. 

The primary objective of our studies was to observe 
the scaling behavior of individual clusters as function of 
time. In certain cases one can relate the scaling exponent 
governing evaporation of small clusters with global scal- 
ing properties of coarsening process due to conservation 
laws [jy} |l2). For example, if the area of small cluster 
S in 2D system evolves as S ~ (to — t) 13 , where to is 
the instant of the cluster's disappearance and (3 is scal- 
ing exponent, and the total number of clusters N evolves 
as N ~ l/t a , for interface-controlled coarsening regime 
the exponents a and (3 are related, a = (3 |Tl] ]. In many 
cases the studies of individual cluster dynamics are more 
reliable than direct studies of global scaling exponents 
which are often masked by the finite size effects. 

The experiments in the mechanical cell were performed 
in the following manner. A uniform gas-like phase was 
prepared by setting the dimensionless acceleration level 
r = Air 2 f 2 Ao/g (f is vibration frequency, Ao is the am- 
plitude of driving A — Aq sin(27r/i), g is gravity accelera- 
tion) approximately at T = 2.28 — 2. 72g. The acceleration 
was then suddenly dropped to the value 1.00 — 1.05 lead- 
ing to spontaneous nucleation of immobile clusters from 
the gas phase (Fig.|l| c,d). 

The area of a selected cluster S was measured as a 
function of time. If the cluster is not the biggest one 
in the system, eventually it begins to shrink. We find 
that the area S evolves according to a power law with an 
exponent that depends on the vibration frequency. For 
frequencies less than 35Hz, S ~ (to — t)P with (3 — 0.63 — 
0.7, see Fig.pl For the frequencies higher than 55Hz, 
S ~ (to — iy with /? w 1. A crossover- type behavior 
observed for an intermediate range of frequencies. In 
this case the area evolves as follows: while the cluster is 
above a certain size, its area decreases linearly with time. 
However, when the cluster decreases below a certain size, 
its area decreases quicker, as S ~ (to — t) 13 with (3 ~ 2/3. 

Qualitatively similar behavior is observed in an elec- 
trostatically driven system (Fig.||). In the electro-cell, 
the uniform gas was prepared by applying an external 
field of 7 kV/cm before "quenching" the gas below E^- 
For the case of a DC electric field applied, as well as for 
an AC field with the frequency less than 60Hz, we find 
that the scaling exponent (3 w 2/3 (actually the range is 




FIG. 1: Shapes of granular clusters. Upper row shows clusters 
in electrostatically driven system, a) / = 20Hz; b) / = 120Hz. 
Bottom row shows clusters in mechanically vibrated system, 
c) / = 25Hz; d) / = 95Hz. Left column: heap-like clusters 
corresponding to 3D behavior. Right column: flat monolayer 
clusters corresponding to 2D regime. 
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FIG. 2: Cluster area vs. time for mechanical system. Circles: 
/ = 25Hz, P = 0.7. Crosses: / = 35Hz. (3 = 1.03 at the 
beginning of evaporation (see inset) and j3 — 0.66 at the end of 
collapse. Triangles: / = 55Hz, f3 = 1.04. Squares: / = 85Hz, 
(3 = 0.7. Inset: S vs t for / = 35Hz in linear scale 

0.63 — 0.7). Accordingly, the cluster area shrinks linearly 
in time for frequencies above 90Hz. In the range from 
60 to 90Hz, there is a crossover behavior observed with 
these two exponents. 

The apparent frequency dependence of the scaling ex- 
ponents (3 for cluster evaporation can be related to the 
transition from two to three dimensional behavior in the 
granular gas. For the case of DC or low frequency AC 
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FIG. 3: Cluster area vs. time for electro-cell. Circles: / = 0, 
13 = 0.7. Squares: / = 10Hz, f3 = 0.67. Diamonds: / = 20Hz, 
(3 = 0.63. Up triangles: / = 50Hz, (3 = 0.68. Crosses: 
/ = 75Hz, f3 = 1.04 at the initial stage of evaporation (see 
inset) and (3 — 0.67 at the end. Down triangles: / = 100Hz, 
f3 = 1.02. Inset: S vs t for / = 75Hz in linear scale 

driving, the particles have enough time to reach the upper 
plate and the granular gas exhibits 3D behavior. Thus, 
gas particles can fly above the cluster and the exchange 
of particles between gas and solid phases occurs at entire 
surface of the cluster. Consequently, the clusters contain 
more than a monolayer of the particles (Fig. |l]a). Mea- 
surements show that in the electrostatic cell with a gap 
equal to 1.5 mm, cluster thicknesses equal to 4 monolay- 
ers are observed in the central part of the cluster. The 
thickness is limited due to an increase of the electric field 
above the cluster as the gap between its surface and up- 
per plate shrinks. When the field value at the cluster 
surface reaches the second threshold value E 2 , particles 
cannot remain immobile, even in the middle of the clus- 
ter, and will rejoin the gas. As the extent of vertical par- 
ticle motion becomes smaller than the particle diameter 
(for high frequencies), the granular gas becomes two di- 
mensional. In this case, the exchange of particles between 
gas and solid takes place only at the outer boundary and 
clusters have a monolayer thickness (Fig [l}d). 

Similar transition occurs in mechanically vibrated cell. 
If we keep the value of the dimensionless acceleration T 
fixed, lowering of the frequency implies shaking at higher 
amplitudes, and, consequently, increase of the extent of 
vertical motion. Thus, at lower frequencies the vertical 
particle jumps will eventually exceed the particles size d, 
i.e. the motion is effectively more 3D than it is at higher 
frequencies. The dimensional transition in the collisions 
dissipation rate has been previously observed in Ref. jl3) . 

The crossover frequency can be determined from the 
comparison of the height of particle jumps £ and the 
particle size d. In mechanical system the height may 
be roughly estimated to be equal to vibration amplitude 
Aq, i.e. £ = g/4ir 2 f 2 since clusterization occurs when 



r«l. From the condition £ « d one find the crossover 
frequency l/2iryj g/d — 40Hz for d — 150/im, which is in 
a good agreement with experiment. The value of £ de- 
creases as the frequency increases and when it becomes 
smaller than a particle diameter, the granular gas ex- 
hibits 2D dynamics. In this case clusters are flat (Fig. 
|l|d). As the frequency decreases, the magnitude of ver- 
tical jumps increases and the clusters assume heap-like 
form (Fig. |l|c). Their evaporation will reflect features 
of 3D gas dynamics as is the case of electrostatic excita- 
tion. The saturation of the cluster occurs when the clus- 
ter hight achieves the magnitude of the vertical jumps 
of isolated particles. It is necessary to mention that this 
heap-like form of the clusters is not necessarily related 
with the interaction between particles and trapped air 
flow through the bulk of the medium as it was in Refs. 
|]l4j| . Clusters keep heap-like form even in the evacuated 
cell (pressure is 5mTorr). 

According to the model presented in the |l5| when the 
cluster size becomes small (i.e. near the moment of evap- 
oration) the normal velocity component of the interface 
between solid and gas phase for 2D dynamics is governed 
by equation R « —K/2R. Here R is the local radius of 
the interface, K is surface tension coefficient which can 
be expressed through the diffusivity of particles. Thus, 
for the cluster area S — ttR 2 one derives S — irK(to — t), 
i.e S decreases linearly during cluster evaporation. 

To match the low-frequency scaling one assumes that 
the evolution of small clusters is governed by equation 
R w —K/R 2 . It yields R 3 ~ (to — t) and the area occu- 
pied by the clusters scales as S ~ (to — t) 2 / 3 . This result 
can be interpreted as a manifestation of the decrease of 
the surface tension coefficient K in multilayered clusters. 
One expects that the coefficient K is inversely propor- 
tional to the cluster's thickness h (it requires more time to 
evaporate a multilayer than a monolayer), i.e. K ~ h . 
Since for heap- like clusters the thickness h is proportional 
to the radius R, one obtains R ~ R~ 2 . This argument 
applies both to mechanical and electo-cells. 

We would like to emphasize that the low-frequency ex- 
ponent j3 — 2/3 appears to not be related to the global 
scaling exponent a for the number of clusters N . This is 
due to saturation of the height of large clusters: when a 
cluster is large enough (R ^ h) the dynamics of its evap- 
oration has a 2D character, as the gas/solid exchange oc- 
curs mostly at the edge of the cluster. However, when the 
cluster size decreases to R ~ h, evaporation evaporation 
occurs on entire surface of the cluster. 

We find that the amplitude of driving mostly affects 
the pre-exponential factor in the scaling law, i.e. the 
surface tension. Some selected results are shown in Fig. 
^. As one sees from the figure, the value of the surface 
tension (and, therefore, diffusivity coefficient for the par- 
ticles) increases with the increase of the applied electric 
field. This behavior can be attributed to the increase of 
the granular gas temperature. 
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FIG. 4: Cluster area vs. time for electro-cell. / = 20Hz. Cir- 
cles: electric field E = 230V7mra. Squares: E = 240V '/mm. 
Diamonds: E = 245V/mm. Triangles: E = 250V/rara. 
Crosses: E — 260V/mm. Inset: the dependence of cluster 
surface tension on the applied electric field. 
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FIG. 5: Current vs. time during coarsening in DC electric 
field. White line shows the signal averaged over 100 points. 
Current peaks correspond to evaporation of clusters. The 
signal/noise ratio is better than 3. 

One of the interesting predictions of the phenomeno- 
logical model is a fine structure in dependence of 
granular gas density vs. time during the coarsening. 
A small cluster shrinks so quickly at its evaporation 
point that the gas concentration temporarily increases 
Ji~5| , pi |i~7[ resulting in a peak feature. The electric cur- 
rent through the cell is carried by the particles belonging 
to the granular gas phase and is, therefore, proportional 
to the amount of gas. This leads to a transient increase 
of the current through the system at the moment of small 
cluster evaporation. These peaks of current were indeed 
observed for the case of electrostatic cell with applied 
DC voltage (Fig. 1st) . Over time, the current approaches 



a value corresponding to the coexistence of one cluster 
and the granular gas. The observation of the peaks rein- 
forces the validity of the continuum model described in 

In conclusion, we studied experimentally the cluster 
dynamics of evaporation in two driven granular systems, 
mechanical and electrical. The transition between 3D 
and 2D behaviors was found. These results are qualita- 
tively similar for both systems. For the cluster coarsening 
process in a granular gas excited by a DC electric field, 
the transient increase of current through the system at 
the moment of a small cluster evaporation was observed. 
We are grateful to B. Meerson for fruitful discussions. 
This research was supported by the US DOE, grant W- 
31-109-ENG-38. 
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